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Abstract: The "Newtonian" or non-relativistic decomposition of Einstein's gravitational 
field is useful in the post-Newtonian approximation. We obtain the full non-quadratic 
Einstein-Hilbert action in terms of these fields as well as the harmonic gauge fixing term 
and find fairly simple expressions. We discuss alternatives to the harmonic gauge. 
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1. Introduction 

Newton's universal law of gravitation can be expressed through the action of a (static) 
gravitational potential (f> 

s = -^g/* a N 2 (L1 » 



m 2 



m0(r 



together with its coupling to massive objects f dt 

In Einstein's theory of gravity (General Relativity, GR) [|| the gravitational field is 
promoted to a space-time metric g^ u subject to the Einstein-Hilbert action 

Seh = --^ [ ^d A xR[g] , (1.2) 



16vrG . 

where the overall sign in ( |1.2[ ) reflects the signature convention for g^ u which we take to be 
time-like. The coupling to massive objects becomes — m f dr , where x^{t) is the particle's 
trajectory, the proper time is defined by dr 2 := dx^ dx u , and for clarity we used c = 1 
units. 

In the Newtonian limit Einstein's gravity reduces to Newton's. Yet, Einstein's theory 
calls for a 10 component field, while Newton's theory has only one. A natural question 
arises: what is the physical role of the other metric components in the Newtonian limit? 

In 2007 we proposed a change of variables from the metric to certain non-relativistic 
gravitational (NRG) fields based on a temporal Kaluza-Klein reduction. One of them 
can be identified with the Newtonian potential </>, and the others also have a role, even 
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if subleading, in this limit, as we review below, and hence they can be called "Newto- 
nian" . However, the action (and hence also the equations of motion) were not fully known. 
Determining the full action including gauge fixing terms is the objective of this paper. 

We carried the computation using a non-orthonormal frame within Cartan's method, 
namely a hybrid method which incorporates both a non-trivial frame and a non-trivial 
metric. The process turned out to be related to existing work in the Kaluza-Klein literature 
|27| , |28| but it was not used so far in the gravitational post-Newtonian context to the best of 
our knowledge. We also computed the harmonic gauge fixing term for these fields, thereby 
completing the evaluation of the total action. 



Background. The Newtonian limit of Einstein's gravity was an issue since its very begin- 
ning |2], |3j, continued with Q and developed into |5[. The interest further arose in the last 
three decades with the development of the Post-Newtonian (PN) approximation, see [||, ^ 
for recent reviews and g | for some earlier work. The central application of the Post- 
Newtonian approximation is the analysis of the two body problem in GR and the resulting 
gravitational wave signature. This problem is crucial for the global effort to directly detect 
gravitational waves [pi] and eventually interpret their signals. 



In 2004 the effective field theory (EFT) approach to GR [11] was put forward. This 
approach borrows from the ideas of effective quantum field theories. It is based on a hierar- 
chy of scales in the problem and essentially replaces the traditional method of dealing with 
finite size objects through matched asymptotic expansion by introducing instead effective 
interactions of point particles with their background. See for early precursors of the 
EFT method. 



In [13, 14] we introduced the NRG fields 1 and used them to give what is probably the 
shortest derivation of the leading post-Newtonian correction, known as the Einstein-Infeld- 
Hoffmann interaction Their utility received strong support from [|l6[ who reproduced 
the 2PN effective action (see || and references therein) through the EFT approach. The 



authors of [16] compared the NRG fields against the standard metric, and found the NRG 
fields to be better suited for the calculational purpose. After the first arXiv version of 
this paper appeared further support was received when the computation was successfully 
extended to reproduce order 3PN using an impressive computerized calculation |p7| . The 
computations in [|l4|, [l(| 17] required the evaluation of certain terms in the NRG action and 



those were evaluated in a perturbative weak field expansion and sometimes with the help 
of computerized computation. The full and concise non-perturbative expressions presented 
here allows to readily expand and read off the gravitational vertices for the PN perturbation 
theory and it should be useful for further research including advancing the post- Newtonian 
state-of-the-art and reaching order ^PiV (which is not available yet in any method), a 
goal announced in [jl?]]. In such higher orders the recursive diagrammatic relations of ]l8| 
should be helpful as well. Related interesting and relatively recent work on PN and/ or 
EFT appeared in flSj| . 

This paper is organized as follows. In section we review the definition of NRG fields. 



1 See |], [l^] for early precursors of the NRG fields. 
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In section || we present the Einstein-Hilbert action in terms of these fields (the derivation 
is relegated to an appendix). In section || we discuss possible gauge choices from the PN 
perspective and compute the harmonic gauge fixing term in NRG fields. 



2. NRG fields 



The non-relativistic limit is defined as the limit of slow velocities (relative to the speed of 
light). As a zeroth approximation we consider the stronger time-independent (stationary) 
limit. In this limit it is natural to perform a temporal Kaluza-Klein (KK) dimensional 



reduction [13, [lj] (the original and standard KK reduction [2C] was spatial, of course 



i\2 



{dt - Ai dx % ) 



-24> 



7ij dx l dx 3 



(2.1) 



This relation defines a change of variables from g^ u to ((f), Ai,jij), i,j = 1,2,3 which we 
call "non-relativistic gravitational fields", or in short "NRG-fields" . 

The action, when translated into NRG-fields and within the stationary limit becomes 



S = - 



1 



16vrG 



dt drxy/j 



-R[-y] + 2 



(2.2) 



V0 



di<j)dj<fi and F 2 is conventionally defined by F 2 = FijF^, F{ 



ij 



where 
diAj - 

The above-defined NRG fields (|2.1| ) resemble the definition of the well-known ADM 
fields (see also the review |22] and references therein) given by ds 2 = N 2 dt 2 - 
1ADM ij(dx l + N l dt)(dx^ + N^dt), where N,N l are the lapse and shift. Moreover, for 
Ai = = N l the definitions coincide as they do at linear order for all fields (after a suit- 
able Weyl-rescaling of the ADM fields). However, there is a marked difference between 
the two: the ADM fields are designed for the initial value problem and the time evolution 
of the metric, while NRG fields are designed for the post-Newtonian evolution of say a 
two-body problem. The NRG fields are defined by a temporal KK reduction, while ADM 
is nothing but a spatial KK reduction (in order to focus on the time evolution), namely 
in NRG the shift is applied to time and in ADM to space. More concretely this means 
that NRG fields transform nicely under time-independent coordinate reparameterizations, 
while ADM transforms nicely in the space-independent case. So in this sense ADM and 
NRG are sort of opposites actually. 

After the first arXiv version of this paper appeared a detailed and practical comparison 
between NRG, ADM and a modified ADM fields was made Q in the PN context and NRG 
was shown to be the more economical field definition starting at order 2PN. 

It is interesting that even though the Kaluza-Klein (KK) ansatz was available since 



1921, it was not applied in full to the post-Newtonian approximation before [13, 14] as far 
as we are aware. One possible reason is that the usual KK procedure is for a compact 
space-like dimension while here it is applied to time. 



Physical meaning of NRG fields 4>,Ai, jij. Comparing the kinetic term for <j> in the 



action (2.2), with the Newtonian field action ( |1.1| ) together with the respective interaction 
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terms we find it natural to identify (f> with the Newtonian potential. To leading order around 
flat-space time (Newtonian limit) this definition coincides with the accepted definitions of 
the Newtonian potential i^tv in the literature, for instance goo ~ 1 + 20jv- Due to the 
nearly stationary nature of the Newtonian limit, we find eft to be a good way to extend this 
definition beyond the linear order. 

The vector potential A has an action which apart for an overall sign resembles the 
magnetic part of the 4d Maxwell action, and accordingly it is natural to call F the gravito- 
magnetic field and call A the gravito-magnetic vector potential. This name originates in a 
certain similarity between gravity and electro-magnetism. The strong similarity between 
Newton's gravitational force and Coulomb's static electrical force, together with the obser- 
vation that the transition from electro-statics to electro-dynamics requires to supplement 
the scalar electric potential by a vector potential, promoted already in the 19th century 
suggestions to add a vector potential to the gravitational degrees of freedom. In fact, it is 
known how to obtain such a vector potential in the weak gravity/ Post-Newtonian approx- 
imation to GR, a point of view known as "Gravito-Electro-Magnetism (GEM)" (see for 
example |24| , |i~5[ ]) and references within). Just like in electromagnetism the vector poten- 
tial couples at leading order to the charged current, which in gravity is the mass current, 
namely S D m J dtv ■ A, and it is responsible for a gravito-magnetic force between two mass 
currents. A spinning object consists of a dipole of mass current and accordingly it couples 
universally to the magnetic field B = V x A through S D ^ J dtJ ■ B. We note the reversed 
sign of the kinetic term for F. This is directly related to the fact that in gravity both the 
current-current and the spin-spin force have an opposite sign relative to electro-dynamics, 
and so "north poles attract" p5|] . 

Finally, starting with order 2PN we must also account for the 3-metric tensor 7^ which 
comes with a standard Einstein-Hilbert action in 3d (this is achieved through the Weyl 
rescaling factor in front of 'jij in the ansatz). This field has no electromagnetic analogue. 

3. NRG action 

In this section we go beyond the stationary approximation of the last section and derive 
the full, time-dependent action for the NRG fields. For greater generality we work in an 
arbitrary space-time dimension d. We perform the change of variables in two steps. In the 
first we replace 





and in the second step 



The first step is achieved through a dimensional reduction 



ds 2 = e 2(t> (dt - Ai dx 1 ) 2 - ^y ij dx i dx j 
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The Einstein-Hilbert action S = 



-l/(16vrG) / RdV becomes 



1 



- je 2 *F 2 - \ e" 2 * (%%n' k T" ~ ^"%)') ~ m) (3-3) 



where we define 



A := di + Ad t 

Fij := DiAj - DjAi = F %j + A t Aj - AjAi , (3.4) 

a dot denotes a time derivative and Rffi] denotes the Ricci scalar of the spatial metric 
jij where the derivatives in its expression are replaced everywhere as follows di — > Di. 
Borrowing notation from the Mathematica software this definition can be stated by 

Rfil : = R[j\ /. di^Di. (3.5) 



Note the following alternative forms for two of the terms in (3.3): lij^lkll^l^ 1 



-l 3 lij an d l lJ lij = 9*logdet7. 

The action (|3.3| ) for the dimensionally reduced fields ( |3.2D is the essential step towards 
incorporating time dependence into the action, and we find that it leaves quite a compact 
result: the F 2 , R terms generalize appropriately the well-known stationary action, while 
the remaining term is a kinetic term for 7^ which defines an associated metric on field 
space, namely on the space of metrics, and it coincides with one of the deWitt metrics p6| . 
In order to compute the action we used a non-orthonormal frame within Cartan's method, 
namely a hybrid method which incorporates both a non-trivial frame and a non-trivial 
metric (see appendix |A|). 

The second step is a Weyl rescaling 



lij = e-^S; (3.6) 



where we use a shorthand notation 



d:=d-3. (3.7) 

Altogether the metric reads now 

ds 2 = e 2 ^ (dt - Ai dx l ) 2 - e-^jijdx^aP (3.8) 



which agrees with ( |2.lD for d = 4. The Weyl rescaling is designed such that the stationary 
action for 7^ will be the standard Einstein-Hilbert action, without the prefactor present 
in (3.3). We denote the perturbation of the spatial metric jij by 



a ij := Hj ~ <% ■ ( 3 - 9 ) 
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Substituting into the action we obtain our main result, the fully non-linear and time- 
dependent Einstein- Hilbert action for the NRG fields 



l (l + V) \Di4>\ 2 - IfWtof* _ R[l] 



+2(1 + l/d)^A i D j 

J — 9 f 14-1 IJ\rh / - . ih. <il I in- \2 

~4 



e 2(1+1/d ^(%-7H7 4fc 7 j7 -(7 iJ ^) 2 ) 
-2(l+l/d)0 ((,, - (^+ 2 )(^+ 1 ) 12 



(1 + -.) 07 u %- " V £ ' Fj I (3-10) 
In particular for d = 4 space-time dimensions the full NRG or Post-Newtonian action 

is 

| 2 |Z^| 2 - I e 4 ^ 2 - %] + 47« 

-J ^ (^WV* " (7 ii 7y) 2 ) " e" 4 ^ (I'o- '^,, - 60 2 ) j . (3.11) 

We tested our expression for the action in several limits. In the stationary limit it 
reduces correctly to (^^)- A strong test is provided by the sector with two time derivatives. 
Since it contains no spatial derivatives, the same sector can be computed by a performing 
dimensional reduction over space, defining fields similar to the ADM decomposition. The 
action is given by an the analogue of (^]^) which is then followed by a change of variables 
to NRG. The result matches with the expansion of the two time derivatives of (|3.3| ) given 
by (H). 

The expression for the action, together with the gauge fixing term (^^) conveniently 
encodes the Feynman rules for the diagrammatic computation of the two-body effective 



action in the EFT method O, M, 16, 18 [. The relevant terms up to order 2PN include 



the quadratic time-dependent vertices 

2 , A 2 (3.12) 
and the non-quadratic vertices of the form 

<jxj) 2 , V(p<j)A, (p(VA) 2 , dijdi^dj^ . (3.13) 



Our expression, together with the recursive relations of [18] should be useful for computing 
higher PN orders. 

3.1 Massive Kaluza-Klein action 

Suppose we take the KK direction to be space-like (as in the original KK idea) and denote 
it by z. Then our action ( |3.10| ) is related to the full action including z-dependent fields, 
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namely the massive modes. It turns out that this action is essentially known in the liter- 
ature: the 4+1 pre-Weyl rescaling action in compact form (analogous to our ( |3.3| ) ) was 
found by Aulakh and Sahdev [27| and an expanded 4+1 post Weyl rescaling expression 
(analogous to our ( |B.1[ ) ) was found by Maheshwari [28] (see also the reviews ]29[]). Below 
we illustrate how to use our result to obtain the full 4 + 1 action after Weyl rescaling and 
in compact form. 

The metric ansatz becomes 

ds 2 = e- 2,f>/ % u dx^dx u - e 2 ^ (dz - dx^ f (3.14) 

where fj,, u = 0, 1, . . . , d — 2. We observe that the substitution 

t O z (3.15) 

amounts to reversing the Ricci scalar (since it reverses the metric in the non-orthonormal 
frame) 

R[ gixv ] ->• -R[g„ v ] • (3.16) 

This allows us to obtain the action by reversing fl3.10p . For concreteness let us present the 
form of the KK action for d = 5 

1 f „ 

S = — / J — 7 d x dt dz 

16ttG J 

!^\D,cf>\ 2 -^F 2 -R[ 7 ] +^D^d z A v (3.17) 

~\ (l w l U(T {d zllxv ){d zlp( ,) - (r u d zl , v f) - e- 3 * Q (^d zl , u )dA - 3(^0) 2 ) 

This expression could be useful for studying non-linear interactions involving massive KK 
modes. 



4. The gauge 

In computing the two-body effective action we are free to choose the gauge (choice of 
coordinates). It affects both the length of the calculation as well as the end result since 
the effective action depends on the coordinate location of the two bodies and as such is not 
gauge invariant. 

The harmonic gauge is used throughout the PN literature. It is defined by 

= := Y>i p g^ = — -L= d u {y/=ggT) (4-1) 

where Ty P is the Christoffel symbol for the space-time metric g. The corresponding gauge 
fixing action is taken to be 

Sgf = J d d x^g^T» (4.2) 
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In the first subsection we evaluate (|4.2| ) in NRG fields (p.8|) . In the second subsection 
we critically discuss the rationale for choosing the harmonic gauge. We note that the 
harmonic gauge is special for being covariant under Lorentz transformations, while the PN 
approximation breaks the symmetry between time and space and therefore choosing the 
harmonic gauge is not obvious. We proceed to analyze this question explicitly order by 
order in the PN expansion up to 2PN. We summarize our results and discuss them at the 
end of that subsection. 

4.1 Harmonic gauge fixing term 



In order to express the harmonic gauge fixing term (4.2) in terms of NRG fields we find in 
appendix [A] how to express in terms of a frame and a frame metric ( |A.16|) . Substituting 
@P we find 



' e P+W 7 « D( ^ +2 (l + 4) e-< 1+1 ^- ie-P+Wg, log 7 Y - |r ( [ 7 ] - A 



where we define 



Fijkh] ■= IW7] /. di -> Di 



= r ijfcM + 2 ( A iiik + Miij - Aij jk ) . (4.4) 

The expression ((O^) was tested by using two additional frame representations of the 
metric: an orthonormal frame and another which shifts the field <j) from the frame to the 
metric. 



4.2 Comparing gauges for the PN limit 

In general, the freedom to choose a gauge could be a blessing provided that it is chosen 
judiciously. Actually given a specific problem there must be an optimal gauge (simply 
because the computational cost should be bounded from below and non-zero). Here we 
are considering the PN approximation, which defines a certain class of problems, meaning 
that only partial information about the system is known. Accordingly we expect that there 
would not be a single optimal gauge, but rather an optimal class of gauges where some 
of the gauge freedom is fixed while some residual gauge freedom remains to fit a given, 
specific problem. 

In this subsection we ask what is the optimal gauge (class) for the post- Newtonian 
approximation. 

The harmonic gauge fixing appears to be standard in the literature. 2 However, the ra- 
tionale for that is not clear to us. One of the appealing features of the harmonic gauge (per- 
haps characterizing it) is its covariance property being a vector under global Lorentz trans- 
formations of (asymptotically) flat space-times. Yet, in the PN approximation, Lorentz 

2 Sometimes one encounters the closely related linearized harmonic gauge. 
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invariance is broken to spatial rotations by the choice of a rest frame (choice of time), and 
hence we proceed to consider a more general class of gauges with less symmetry order by 
order in the PN expansion. 

The Newtonian order. The newtonian field is gauge-invariant to leading, stationary 
gauge transformations. Accordingly the question of gauge choice is moot. 

Quadratic stationary gauge fixing starting from 1PN. In the PN approximation 
we decompose the gravitational field into the Newtonian potential (f>, the gravito-magnetic 
vector A{ and perturbations a%j of the spatial metric tensor 7^ ( p.8|j3.9| ). The vector 
appears in the stationary, unperturbed action (2.2) in a Maxwell form, and accordingly the 
standard Lorentz gauge fixing term must be chosen at the quadratic level 

C = V • A . (4.5) 

in order to avoid the propagator from mixing the different vector components. This gauge 
fixing term affects order 1PN through the diagram of gravito-magnetic interaction (vector 
exchange) . 

Similarly 7«, the spatial metric has an Einstein-Hilbert action and for simplicity of 
the propagator we better choose the harmonic gauge (note: only in the spatial directions), 
at least linearly. This gauge choice shows up first at order 2PN through the diagram of a 
exchange. 

Quadratic but non-stationary sector starting from 1PN. Before gauge fixing 
the action includes a term of the form A ■ Vc/> just like the post-Coulomb approximation 
to electro-magnetism. This term induces a mixing 2-vertex for eft, A. A generic gauge 
fixing term includes also a <fi 2 vertex and hence generically three diagrams contribute: an 
exchange of eft with a retardation vertex, a <f> — A exchange diagram and a diagram with a 
<j> — A — (j> chain. 

However, after integration by parts the <j> — A 2-vertex is proportional to (j)V 'A and 
can be removed by shifting the gauge condition on V^4, more specifically by deforming the 
gauge condition (|4.5|) to become 

C = V-1+ 2^1 + 1^0. (4.6) 

This constraint turns out to be part of the harmonic gauge and hence up to order 1PN the 
harmonic gauge is certainly optimal. We note that just like the electro-magnetic two-body 
effective action is gauge invariant at any post-Coulomb order, the same holds in the first 
order of the gravitational case, due to the similarity of the actions up to this order, and 
hence the sum of the three diagrams is independent of gauge. 

Similarly, at order 2PN with harmonic gauge the diagram of A exchange with a retar- 
dation vertex replaces two diagrams with a a A mixing which would appear in a generic 
gauge. 

Non-Quadratic gauge fixing starting from order 2PN. We still have freedom 
to deform the gauge constraint fl4,6|) with non-linear terms. This freedom can be used to 
cancel any term in the action which is proportional to the linearized gauge constraint. 
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Figure 1: Diagrammatic representation of the new static interaction obtained after gauging away 
A l Wi66 from the Einstcin-Hilbert action. 



P 

► — « 



k, i y 



k + P,J ^ s 




(b) 



Figure 2: New static Feynman diagrams obtained after gauging away A l Vi4>(j) from the Einstein- 
Hilbcrt action, (a) and (b) replace figures 5(1) and 5(k) of Gilmore-Ross respectively. 



At order 2PN the vertex VcfxfiA appears in two diagrams: 5(m) and 5(n) of 16]. We 
may cancel this term by deforming the gauge condition as follows 



C n = V-l+2( 1 + i 



+ A l di> 



(4.7) 



At the same time this gauge fixing produces a term of the form V(f)AVA which at first 
sight is already present (see eq. (35) in [fl6|] ), however the tensor structure of the new term 
is different and one must add the static Y-type interaction presented in figure |l[ 
Computing the value of the terms added to the diagrams we find 



FigIa) = ^|(16vrGm 2 ) 2 ^ 



dt 



d d 1 k d d l p (v 2 • p)(v2 ■ k) 



rd-l 



Figi(6) 



d-2 y ' 2 



dt 



(2^-! (2 7 r) d - 1 p 2 k 2 (p + k) 2 

d d ~ l k d d ~ 1 p (vi-p)(v 2 -p) ik ., 
{2Tt) d - 1 (2n) d - 1 p 2 k 2 (p + k) 2 6 



(4.8) 



Integrating over one of the wave numbers and evaluating the resulting Fourier transform 
by means of (A. 2) and (A. 10) of |18[ respectively, yields (in d=4) 



Figi 
Fig.| 



G m^mi 



2(v 2 ■ f) 



3(£i • f)(v 2 ■ r) - (vi ■ v 2 ) 



(4.9) 



We conclude that while in new gauge (4/7) avoids computing two diagrams (fig. 5(m,n) 
of [|l(|) it generates two new terms in the diagrams in fig. 5(k,l) which appear to be of 
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comparable computational difficulty. Hence it is not clear that the new gauge d4.7| ) is any 
better than the harmonic. 



Comments about 3PN. At order 3PN log terms appear, reflecting an RG flow H, 11 



17]. The associated divergences are canceled by available redundant terms on the world- 



line which reflect the choice of coordinates, but not in observables. Presumably this is a 
price for working with harmonic coordinates. Indeed, the equation for changing the radial 
coordinate of the Schwarzschild solution from standard to harmonic has log solutions even 
in asymptotically flat space, and those terms could be essential in background which are 
not asymptotically flat. 



Optimal Post-Newtonian gauge: summary and discussion. We critically analyzed 
the optimal choice of gauge for the PN approximation. We found that the harmonic gauge 
is certainly optimal at the level of propagators and 2- vertices (starting with 1PN). At 
2PN we found a different gauge which reduces the number of diagrams. However, this 
does not lead to a reduction in computational difficulty. The reason for that is that each 
standard Feynman diagram may represent several computational terms corresponding to 
the positioning of derivatives and to various index structures. A faithful diagrammatic 
representation of the computational difficulty requires to distinguish vertices of a field 
from those of its derivatives and also to reflect the index structure. Hence altogether we 
do not find a gauge which is better than harmonic, at least up to 2PN. 

Our analysis suggests that the harmonic gauge could be the optimal PN gauge (in 
the bulk). In hindsight, the reason for that could be simple. The total PN action is 
composed of the bulk gravitational action and the interaction with the object's world- 
line. The computational difficulty is dominated by bulk irreducible loops, which depend 
on bulk vertices while non-linear world-line vertices contribute only factorizable diagrams 
(see for example |]18[). Hence the choice of gauge depends only on the bulk action which 
enjoys an increased symmetry (relative to PN). The harmonic gauge can be characterized 
by its background independence, namely in any background the quadratic action for the 
fluctuations is given simply by the Lichnerowicz operator. Probably this characterization 
leads to its advantage also for non-quadratic orders around flat space. 
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A. Deriving the NRG action 



In this section we shall derive ( |3.3| ) - the full Einstein-Hilbert action in terms of the dimen- 
sionally reduced fields Q3.2| ). We use a non-orthonormal frame within Cartan's method, 
namely a hybrid method which incorporates both a non-trivial frame and a non-trivial 
metric, see for example 3C]. 

The dimensionally reduced fields (|3.2[) define a natural frame of 1- forms 



0* = (dt - Ai dx l ) 

9 l = dx l (A.l) 

and a metric 

g ab = diag(l,-7ij) (A.2) 

where a,b = (t,i) are frame indices and from hereon we shall replace 7 — > 7 throughout 
the derivation for clarity of notation and without risking confusion. 



It is convenient to record the inversion of (A.l 



dt = e~ V + AiP 
dx i = 6 l (A.3) 



as well as the dual basis of vectors 

ie t = e-^d t id t = e^e t 

\e. = d i + Aid t = Di \di = e x - t^Api 

The connection 1- forms uj a b are the solution of 

d6 a + co a b 9 b = 
u a b + u ha = dg ab . (A.4) 

In order to solve this system one defines the components of u a b to be 

u a b = u a hc 6 c (A.5) 
and these components are given by the sum 

Wafec = Sabc + Tafec (A-6) 

where Ti abc is related to the torsion of the frame, T abc is the Christoffel symbol for the frame 
metric g ab (it is symmetric in 6c, namely r a ^ c ] = 0) and both are defined and computed 
below. 

First one defines a tensor C ab by 

d6 c + ^C ab c a 8 b = . (A.7) 
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C ab c is antisymmetric in ab, namely Cy m = 0. In our case one has 



cJ 



e Fij 
Did + Ai 



where F,Di are defined in (3.4). 
S a ft c is defined by 



(c, 



abc 



G 



bca 



a 



cab) 



(Ai 



(A.9) 



S a b c is also antisymmetric in ab, namely £( a ft) c = 0, and in our case we have 



Hit 



% = 2 
Did) + Ai 



'j 



Next one computes the Christoffel symbol defined as usual by 



r a fe c = ^ (db§ac + d c g a b — d a gbc 



(A.10) 



(All) 



In our case it is given by 



r 



jk 



ij 



tj 



(A.12) 



where V is defined by (4^). Substituting back ( [A.10 , A.12 ) into ( A.6 ) provides us with the 
connection u) a b c - 

The Einstein- Hilbert action is expressed in terms of the curvature, which is related to 
derivatives of the connection. However, after integration by parts we found that it could 
be expressed in terms of the connection only, so there is no need to compute the curvature 



(-WttG)S eh = J RdV~ J ^9 (u abc Lo abc -u) a Lu ! 
where ~ denotes equality up to boundary terms and we define 



(A.13) 



u a := g bc u abc = g bc T 



" ba 



abc 



ra I p a 
ba ~r u afe 



Ca 
ab 



(A.14) 



We note that if one were to consider J fdVR where / is any scalar function, the integration 
by parts would add the term — (uo a — td a )d a log /. 

Substituting the values of co a bc into ( [A. 13 ), we encounter the expression 



f^f kij - 1 kl Y i kl + ( 7 «T % - 7 fei f L)((A0 + Ai) 



(A.15) 
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where everywhere the connection here is of the spatial metric alone = r^]^. Recalling 
that the measure is d d ~ l x dt and using the comment at the end of the last paragraph 

this expression can be integrated back by parts to yield —^[7]. Finally replacing back from 
the temporary notation of this section 7 — > 7 we obtain the expression for the action (|3.3| ) 
that we sought to derive. 

Harmonic gauge fixing term. In order to calculate this term it is useful to express 
the harmonic constraint in terms of a hybrid frame formalism. We obtain the following 
relation 

T a [g up ] := e£ T p [g up ] = T a [g bc ] + 1^ (A. 16) 

where the first term is the usual expression only in terms of the metric in the frame basis 

Talgbc] ■= 9 bc (d b g ac - \d a9b A (A.17) 



2 

while the second term involves also the frame 

K := {g ab g up + « - < <) d v 6 b p (A.18) 

B. NRG action expanded 

In this section we present expressions for the expanded action and group it according to the 
number of time derivatives for the purposes of the PN approximation. These expressions 
were computed using a rather different method (relying on a presentation of the metric 
by an orthonormal frame) and were confirmed to be identical with the expanded form of 
(|3l0|) . 



Since each term in the scalar curvature contains precisely two derivatives of the metric 
with respect to either a spatial or a temporal coordinate, one can organize all the terms 
into three groups according to the number of temporal derivatives 

Seh{9) = T^G J R i9]^S^( 1 ,A,cP) + S^( lj A,cP) + S^( 1 ,A, ( j ) ) . (B.l) 



We find 



16vrG 



S (0) (7, A0) = -7^ / d*Xyfy -Rh] + d ~A - 



d-3 K 4 



(B.2) 



1 j 

7 

-2^vV(i* + M 



where a dot denotes the derivative with respect to time, the metric 7 is being used to 
raise and lower the indices such as (<9</>) 2 = 7^ d^djcf) and we use the standard definitions 
F 2 = F {j and i™, Fy = diAj - djA,,. Finally (7, A, (f)) is given by 
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+ 



, d x-Jj 



2" 



Y J AiAjA 2 - ( ^4Mj 



2(d-2)<£/(d-3) 



1 d-2 



+ 
+ 



16vrG d 
1 d 


-3. 


16vrG d 


-3. 


1 / 


d d x 


16vrG J 


1 d- 




8ttG d- 


-,! 



d d x^j ft 



A 2 + d~ 1 g _2(d-2)*/(d-3) 

d — 3 



d d x^e- 2(d - 2W( ^ 3) 07 iJ ^ 



J d d x^/ 



> A l Ai 



(B.3) 



with A 2 = Y j AiAj 
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